Structure and Vibrations of the Vicinal Copper (211) Surface 
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We report a first-principles theoretical study of the surface relaxation and lattice dynamics of the 
Cu (211) surface using the plane-wave pseuodopotential method. We find large atomic relaxations 
for the first several atomic layers near the step edges on this surface, and a substantial step-induced 
renormalization of the surface harmonic force constants. We use the results to study the harmonic 
fluctuations around the equilibrium structure and identify three new step-derived features in the 
zone-center vibrational spectrum. Comparisons of these results with previous theoretical work and 
with experimental studies using inelastic He atom scattering are reported. 
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I. INTRODUCTION 

At a vicinal surface the average surface normal is 
tipped away from a high symmetry low-index direction, 
and the microscopic structure consists of an ordered array 
of atomic steps which separate low-index terraces. The 
microscopic atomic structure near a step edge is very dif- 
ferent from that found in the bulk or at a low-index sur- 
face, and as a consequence one expects a substantial re- 
distribution of the valence electronic charge density near 
the step edges. This can play an important role in mod- 
ifying the local dynamical properties and the chemical 
reactivity. Therefore stepped surfaces are important sys- 
tems to study exprimentally and theoretically. In fact 
the influence of steps on the surface chemistry can be so 
large that the macroscopic reactivity of many real sur- 
faces is governed by dynamical phenomena at steps. Un- 
derstanding these processes at the microscopic level can 
lead to control of surface dynamical phenomena through 
control and design of the surface structure. 

Due to the importance of these structures, the dy- 
namics of stepped surfaces have been the subject of sev- 
eral experimental investigations using modern surface- 
science techniques within the last several years. Witte 
and coworkers have studied the effects of surface steps on 
the surface vibrational spectra of Cu(211) and Cu(511) 
using inelastic He atom scattering ^ . These studies were 
able to resolve clearly the zone-folding effects in the low 
frequency surface vibrational spectra due to the ordering 
of single atom steps on these surfaces. However, this work 
did not address the renormalization of the interatomic 
interactions which occur near the step edges, and the 
authors concluded that their observations could be inter- 
preted within a model in which the parameters describing 
the surface atomic interactions were assigned their bulk 
values. More recently Niu et al. ||^ studied the dynamics 
of Ni(977) using inelastic He atom scattering. This work 
revealed a new step-derived branch in the surface vibra- 
tional spectrum which the authors associated with a new 
kind of surface sound wave, bound to the surface along 



the surface normal, and localized laterally to the sur- 
face step-edge. To understand the dynamics of this step- 
induced branch, these authors introduced a phenomeno- 
logical model in which the step edge modes were inter- 
preted as the modes of a stretched string running parallel 
to the surface. Using this model, the authors found that 
the effective elastic constants describing the harmonic 
fluctuations of the string were strongly renormalized from 
the values appropriate for bulk Ni or for a low-index Ni 
surface. However, subsequent theoretical work, based on 
a continuum model for a semi-infinite elastic medium, 
demonstrated that the step-derived features can be un- 
derstood using continuum dynamics and enforcing the 
surface boundary condition appropriate to a weakly cor- 
rugated stepped surface 1^ . Thus the changes of the sur- 
face interatomic interactions near the step edges remains 
an important but unresolved issue. 

In this paper we report a study of the surface struc- 
ture and dynamics of the Cu(211) surface using first- 
principles pseudopotential theory. Previously, the equi- 
librium structures for various stepped vicinal surfaces 
have been examined theoretically using empirical inter- 
action potentials tight binding theory [Q, effective 
medium theory Q and first-principles density functional 
theory The former three methods are only tested in 
a highly coordinated bulk environment; the transferabil- 
ity of these approaches to the structures found on vicinal 
surfaces is thus open to question. First-principles total- 
energy calculations provide a parameter free method 
for studying equilibrium structures. Unfortunately the 
structural complexity of vicinal surfaces has made them 
nearly inaccessible by first-principles methods, particu- 
larly for the transition metals. Nonetheless several years 
ago Nelson and Feibelman reported the first ab initio 
study of the equilibrium structure of a simple metallic 
vicinal surface, Al(331). These authors found that the 
atomic relaxations near the step edges are significantly 
enhanced with respect to those found at close-packed 
low-index surfaces, although the dynamical properties of 
the Al(331) surface were not investigated in their work. 
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Recent improvements in computational algorithms and 
hardware now allow the detailed exploration of realis- 
tic models for vicinal surfaces of transition metals using 
first-principles total-energy methods. 

In this paper we apply the first-principles pseudopo- 
tential theory to study the equilibrium structure and dy- 
namics of the Cu(211) surface. We find that atomic re- 
laxations for Cu(211) are strongly enhanced in the first 
several layers near the surface step edges. This relaxation 
strongly changes the effective surface harmonic force con- 
stants from the values appropriate for the bulk or for a 
low-index surface. Nevertheless, we find that the dom- 
inant effect on the surface vibrational spectra is simply 
the lowering of the surface translational symmetry and 
the folding of the surface Brillouin zone at the stepped 
surface. Specifically, in our calculations we find three pri- 
mary step-induced features in the vibrational spectrum 
at the surface Brillouin zone center. Two of these are the 
"backfolded" versions of the elementary surface elastic 
waves of the close packed (111) surface, while one is due 
to the renormalization of the surface interatomic force 
constants. One of our calculated low-energy modes is in 
good agreement with an excitation observed by inelastic 
He scattering experiments on Cu(211) |^], whereas the 
symmetry and energies of the other two features made 
them inaccessible in these measurements. In this pa- 
per we suggest other experimental probes which could 
be used to examine these latter excitations. 

The remainder of this paper is organized as follows. 
In Section II we summarize the symmetry and struc- 
tural features of the Cu(211) surface, and in Section III 
we briefly present our computational methodology. The 
equilibrium structure predicted in our calculations is pre- 
sented in Section IV. In Section V we apply the model to 
determine the harmonic interlayer force constants near 
the surface, and use them to examine the surface vibra- 
tional excitations. In Section VI we analyze the step- 
induced features in the surface vibrational spectra. A 
discussion and summary of the results is presented in 
section VII. 



II. STRUCTURE OF THE COPPER (211) 
SURFACE 

The ideal lattice termination of the (211) surface of an 
fee crystal is shown in the top panel of Figure 1. A more 
complete geometric classification of the atomic structure 
for stepped surfaces of fee crystals with close-packed step 
terraces is given in Reference 0. At the (211) surface 
the surface normal is rotated by 19.5° away from the 
(111) direction with the axis of rotation along the (Oil) 
direction exposing (111) terraces containing 3 atoms sep- 
arated by single atom (100) facets at the step edges. The 
Cu(211) surface has a primitive rectangular unit cell cor- 
responding to a 3 X 1 reconstruction of the bulk lattice 
structure. The surface retains a mirror plane symme- 
try perpendicular to the (Oil) direction. Although this 



surface cell exhibits a microscopic faceted structure at 
the atomic scale, the unit cell is relatively compact and 
can be studied theoretically using conventional thin film 
supercell methods. The calculations reported below are 
carried out on a structural model in which a thin film 
containing 17 atoms in the primitive cell is terminated on 
its (211) and (211) planes. An illustration of this struc- 
tural model, with the inequivalent lattice sites identified, 
is shown in the bottom panel of Figure 1. 

III. COMPUTATIONAL METHODS 

The calculations are carried out using density func- 
tional theory within the local density approximation 
(non-relativistic) ||]. The electron-ion potential for Cu is 
modeled using an optimized non-local soft core pseudopo- 
tential for the Cu 11-1- ion. Cu, like other first row transi- 
tion metals, is challening for pseudopotential theory be- 
cause it possesses an unscreened 3d valence shell. This 
implies that the peak in the Sd-derived valence charge 
density is relatively tightly bound to the core, and thus 
a well converged calculation requires a large energy cut- 
off in the plane wave basis set. This difficulty has been 
addressed and solved by the introduction of optimized 
pseudopotentials, which allow a well converged calcula- 
tion of the 3d derived levels with only a modest energy 
cutoff in the plane wave basis set. A discussion of the 
construction and testing of these potentials can be found 
in Reference Q. The reference state for our potential 
is the +1 ionic configuration for Cu with orbital occu- 
pations 4s"'^^4p°-^^3(i^. The one-electron wave functions 
are expanded in a plane- wave basis set, up to a kinetic 
energy cutoff of 50 Ry. Using this cutoff we find that 
the truncation error in the computed total energy is less 
than 1 mRy per atom. In the 17 layer thin film shown 
in Figure 1, the thin film structure is repeated periodi- 
cally along the x direction with a 14.5 A vacuum layer 
separating the neighboring films. For this structure ap- 
proximately 25,000 plane waves are required in the basis 
set for the electronic degrees of freedom. We minimize 
the energy with respect to the plane-wave expansion co- 
efficients using a pre-conditioned conjugate gradient al- 
gorithm described in detail in Reference 

Using this method for the bulk fee structure in cop- 
per we calculate an equilibrium nearest-neighbor sepa- 
ration of 2.554 A, which agrees quite well with the ex- 
perimental value of 2.553 A. As a test of the lattice 
dynamical properties of Cu using this potential, we also 
carried out frozen-phonon calculations of the (100) zone 
boundary transverse and longitudinal acoustic phonons 
in the bulk; the calculated (experimental) energies are: 
TA: 32.2 (30.2) meV, and LA: 22.8 (20.7) meV The 
bulk modulus is 1.63 x 10^^ N/m^ which is larger than, 
although comparable to, the observed value 1.42 x 10^^ 
N/m^ |1^. All these results are in good agreement with 
experiment and in excellent agreement with several pre- 
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vious theoretical studies of the bulk elastic proper ties of 
Cu within the local density approximation |12|-fl^]. Our 
calculations are nonrelativistic, that is we have not in- 
cluded scalar relativistic corrections in our theory. These 
are known to introduce a slight (1 % - 1.5 %) additional 
contraction in the predicted lattice constant, and a cor- 
responding increase in the bulk modulus. 

To study the harmonic fluctuations around the equi- 
librium structure, we compute the distribution of forces 
at all of the lattice sites of the thin film produced in 
response to separate small displacements of each of the 
atoms away from their equilibrium locations. The maxi- 
mum ionic displacement that we studied in these calcula- 
tions is 0.05 A. Using the Hellmann-Feynman theorem, 
the induced force on the i-th site in the a-th polariza- 
tion is Fi^a = —dU/dxi^a — —{dH/dxi^a) where H is 
the Kohn-Sham Hamiltonian for the coupled electron-ion 
system, and the brackets denote an expectation value in 
the relaxed electronic ground state in the presence of the 
ionic displacement. The harmonic interaction constant 
coupling the ia and j/S displacements is Ka^p{i,j) = 
d'^U /dxiadxjp = —dFia/dxjp . The displacement fields 
and forces calculated by this method describe the dy- 
namics at the center of the folded Brillouin zone for the 
stepped surface, and so in the analysis below we consider 
only the vibrational spectra at T . 

The harmonic interlayer force constants calculated 
from our thin film describe the physics only within the 
region near the surface (half the thickness of our film). 
This is insufficient to describe realistically the dynamics 
of a semi-infinite crystal terminated on the (211) plane. 
To remedy this problem, we carry out a separate calcu- 
lation to obtain the analogous interlayer force constants 
for Cu in its bulk fee structure. We then construct a 
composite model in which the surface region is recoupled 
to a very thick unperturbed bulk film containing « 1000 
atomic layers. This expanded film is sufficiently large 
to allow us to locate critical points in the bulk phonon 
spectra accurately and to represent the bulk and surface 
modes of our system properly. Since the dynamical ma- 
trix is computed at the center of the surface Brillouin 
zone, mirror reflection through the (Oil) plane remains a 
valid symmetry, and therefore the dynamical matrix can 
be decoupled into two sectors describing lattice modes 
which are even and odd under reflection through this 
mirror plane. 

IV. EQUILIBRIUM STRUCTURE 

Figure 2 illustrates the atomic structure of the (211) 
surface of our fee crystal. The arrows represent magnified 
displacements of each of the atoms from its position in the 
ideal termination of the bulk Cu lattice to its position in 
the equilbrium structure (for clarity these displacements 
have been magnified by roughly a factor of 20) . By sym- 
metry the atomic displacements in the relaxed structure 



are confined to the plane perpendicular to the step edge 
(i.e. the (Oil) plane). We find that the lateral displace- 
ments, i.e. along the (111) direction, while allowed by 
symmetry, are all relatively small; the largest of these 
corresponds to « 1.5% of the lateral spacing in the bulk. 
By contrast the displacements along the surface normal 
are quite substantial near the surface. The relaxation of 
the distance between the two outermost atomic layers is 
— 14.4% of the bulk interlayer spacing in this direction. 
For the fully relaxed geometry the fractional change of 
the nearest-neighbor distance from the step edge atoms is 
nearly 4% of the bulk nearest-neighbor spacing. For com- 
parison, at the Cu(lll) surface the inward relaxation of 
the top surface layer reduces the nearest-neighbor spac- 
ing by 0.8% and for Cu(lOO) the nearest-neighbor spac- 
ing is reduced by 1%. Thus for the vicinal surface the 
atomic relaxations along the surface normal are substan- 
tially larger than obtained for a close packed low index 
surface. The relaxation within the surface plane is quite 
small. 

The displacement pattern shown in Figure 2 illustrates 
the Smoluchowski effect which predicts a reduction 
of the height of the step edge in the equilibrium structure. 
Qualitatively the conduction electron density leaks away 
from the crystal at the step edge, tending to smooth the 
step edge profile. It is easy to see that this generates an 
electric field which pulls the ion at the inside corner of 
the step edge away from the thin film and pushes the ion 
at the outside corner towards the film. A detailed tab- 
ulation of the calculated atomic displacements is given 
in Table L In addition these results show that the struc- 
ture obtained near the center of the slab is essentially the 
bulk fee structure, and so the film thickness is sufficient 
to describe the relaxation of the isolated (211) surface 
accurately. 

V. INTERLAYER FORCE CONSTANTS 

To study the harmonic fluctuations around the equilib- 
rium relaxed structure, we calculate the interlayer force 
constants in our film following the method outlined in 
Section III. The atomic relaxation at the surface leads to 
very large modifications of the harmonic interlayer force 
constants near the surface region. In Table II we report a 
set of selected interatomic interactions near the surface, 
and compare these with their corresponding unperturbed 
bulk values. We can clearly identify the following trends. 
The interlayer coupling constants for displacements po- 
larized along the surface normal is much larger at the 
stepped surface than in the bulk. For example the har- 
monic interaction constant which couples the displace- 
ment of surface site 1 to the nearest-neighbor subsurface 
site, Kxx{^,'^)^ is 4.52 eV/A^ at the surface while the 
comparable interaction is 2.89 eV/A^ in the bulk. The 
intra-terrace interactions are also strongly renormalized 
at the surface; as an example we find that Kxx{^,'^) = 
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0.41 eV/A^ at the surface terrace, which is to be com- 
pared with the bulk value 0.12 eV/A^. For displacements 
polarized within the surface plane, the surface force con- 
stants are again larger, although the changes of the force 
constants are more modest, and quite comparable to 
those obtained for a flat low-index surface |17|. For ex- 
ample the surface (bulk) interaction constant Kzzi^,"^) 
= 3.24(2.70) eV/A2 and if^^(l,4) = -0.49(-0.21) eV/A^. 
Note that the latter represents a very large fractional 
renormalization of the surface interaction; this is mainly 
due to the fact that the interlayer coupling between dis- 
placements parallel to the surface plane is relatively small 
in both environments. The renormalization of the sur- 
face interaction constants for displacements parallel to 
the step direction is similar, e.g. Kyy{l,2) — 1.27(0.92) 
eV/A2 and Kyy{l,A) = 1.47(0.92) eV/A^. Generally we 
find that all the interaction constants are strongly en- 
hanced near the step edges at the surface and this effect 
tends to offset the rather large decrease in the coordina- 
tion number for the surface species. This is most clearly 
seen by inspecting the diagonal interaction constants for 
atoms on a common terrace. For example, the interac- 
tions Kxx{iii) for the three sites i at the surface terrace 
have the values -5.77, -4.61 and -4.37 eV/A^ in the re- 
laxed structure for sites i=l,2,3. Site 1 is the site at 
the outside step corner with the smallest coordination, 
and site 3 is the site at the inside corner with the largest 
coordination for a terrace site. 



VI. SURFACE VIBRATIONAL SPECTRA 

We combine the interlayer force constants for the 
Cu(211) surface region with the values calculated for bulk 
Cu to obtain a composite dynamical matrix for a thick 
Cu film which terminates on a (211) lattice plane. The 
dynamical matrix is decoupled into two sectors describ- 
ing the modes of the combined systems which are even 
and odd, respectively, under reflection through the (Oil) 
mirror plane. For a slab containing N atoms, the di- 
mensions of the even and odd subspaces are 2N and N 
respectively. 

In Figures 3 and 4 we present results for the symmetry- 
projected densities of vibrational states as a function of 
frequency at the F-point in the surface Brillouin zone. 
Figure 3 gives the density of modes which are even under 
reflection through the surface mirror plane. The thin 
solid curve is obtained from a sum over all the modes of 
the composite slab, while the bold solid curve gives the 
ZocaZ density of states ns{uj) = '}2nl2i,a l0ia,nP<5(t^-Wn) 
where 4>ia^n is the amplitude of the n-th normal mode 
on site i with polarization a, and where the sum over 
ia runs over the symmetry allowed displacements of the 
three atoms on the outermost surface terrace, and the 
sum over n runs over all the modes of the slab. 

The surface projected density of even states clearly re- 
veals two step-induced features in the surface vibrational 



spectrum at f . The lower of these, at an energy of 11.5 
meV, is a backfolded image of the Rayleigh wave of the 
"flat" (211) surface which would appear at a momentum 
q = 1.005 A-i along the (111) direction in a continuum 
theory We flnd that the g = ± 1.005 A^^ backfolded 
branches are strongly mixed by the cell periodic varia- 
tions in the surface dynamical matrix so that only the 
lower branch remains localized below the nearby bulk 
transverse threshold at 12.4 meV. This excitation is a 
surface "standing wave" trapped on the surface terrace. 
This can be seen more clearly in the map of the displace- 
ment pattern for this feature plotted in Figure 5. It shows 
that the largest atomic displacements occur at the center 
of the surface terraces and are reduced at the step edges, 
an effect which can be seen clearly in the subsurface dis- 
placement held. Nevertheless this mode possesses a large 
amplitude at the surface, polarized along the surface nor- 
mal, and so this branch is expected to be strongly allowed 
in inelastic atom scattering from the surface. Indeed, the 
calculated frequency agrees closely with the 10.9 meV 
feature observed in inelastic He scattering experiments 
from Cu(211) |l|l. 

The backfolded Rayleigh branch overlaps the projected 
bulk continuum at F, yet the mode remains quite sharp 
in the surface-projected vibrational spectrum. Within 
a continuum elastic theory, mixing between the back- 
folded Rayleigh mode and the projected bulk modes is 
dynamically forbidden at F, since in the continuum the- 
ory the coupling between the two is mediated by the 
surface strain which vanishes for the lowest branch of 
bulk acoustic waves at F This dynamical symme- 
try can be weakly violated in a lattice model. However, 
the sharpness of this surface-derived feature in the vi- 
brational spectrum shows that the mixing remains very 
weak, even for the relatively narrow terraces found on Cu 
(211). ^ 

At higher energy we flnd a second strong step in- 
duced surface feature at 24.7 meV which corresponds 
to a backfolded image of the surface longitudinal reso- 
nance of a flat surface. The displacement held associated 
with this mode is illustrated in Figure 6. Note that the 
surface displacements have their largest amplitudes po- 
larized mainly parallel to the plane containing the sur- 
face terraces. This mode can also be understood as a 
surface standing wave which is trapped laterally by re- 
peated scattering from the step edges. However, it is ap- 
parent from the plot of the displacement field, that the 
effect of scattering from the step edges is much weaker 
than one finds for the lower energy backfolded Rayleigh 
mode. This is because this higher energy feature is not 
a true surface-localized mode, but is resonantly mixed 
with propagating bulk acoustic waves which penetrate 
infinitely deeply into the bulk. Since the displacement 
field for this mode is polarized mainly within the sur- 
face plane we expect it to be a poor candidate for study 
by inelastic atom scattering, which is most sensitive to 
atomic displacements along the surface normal. Never- 
theless, it is very interesting that, unlike the situation 
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for a flat low index surface, a linear coupling between 
the backfolded longitudinal resonance and charge fluc- 
tuations along the surface normal is allowed by symme- 
try in this system since both excitations are even with 
respect to the (Oil) mirror plane. This opens the pos- 
sibility of probing this excitation by inelastic scattering 
of a charged probe (e.g. high resolution inelastic elec- 
tron scattering) from the surface. Further theoretical 
work is needed to quantify the magnitude of the charge 
fluctuations coupled to this mode. We expect that the 
important charge fluctuations coupled to the zone folded 
longitudinal resonance will be found at the step edges, 
and that this coupling can be used to probe the dynam- 
ics of the conduction electrons near the steps. 

Figure 4 shows the zone center vibrational spectra 
for modes that are odd with respect to the (Oil) mir- 
ror plane. The thin solid curve represents the spectrum 
summed over all the modes of the film, while the bold 
curve gives the modes weighted by the squares of the 
displacement fields in the surface terrace. Here the step 
induced features, though present, are much more subtle. 
A broad transverse resonance, labeled TA2, is found cen- 
tered at 11.3 meV just below the bulk projected critical 
point. The displacement field for this excitation is plot- 
ted in Figure 7. We see that the mode corresponds to 
a transverse shear wave which is strongly laterally con- 
fined at the step edges. This mode is completely absent 
within a simple continuum elastic theory, and its appear- 
ance in the spectrum directly reflects the rather substan- 
tial renormalization of the interatomic force constants 
within the surface region. This mode is a poor candidate 
for study using neutral-atom scattering since, as shown 
in Figure 7, it is polarized in the surface plane, along 
the (Oil) direction and is odd under the surface mirror 
plane. Likewise linear coupling between this mode and 
z-polarized charge fluctuations is symmetry forbidden, 
and therefore this mode will be silent in charged-particle 
scattering as well. 

VII. DISCUSSION 

In this section we compare our results to experimen- 
tal measurements of the surface phonon spectra on the 
Cu(211) surface and to related previous theoretical work. 

Witte et al. |l] used inelastic He scattering to study 
the phonons of the Cu(211) surface. These authors ex- 
perimentally studied the surface vibrational excitations 
below 15 meV as a function of momentum transfer in the 
surface plane, both parallel and perpendicular to the sur- 
face steps. At the surface Brillouin zone center, F, an ex- 
citation polarized along the surface normal at an energy 
of 10.9 meV was observed and assigned to the backfolded 
image of the Rayleigh wave of a flat surface. This agrees 
quite well with the Rayleigh mode we calculate at 11.5 
meV. The experiments do not extend to sufficiently high 
energy to detect the backfolded LA resonance, which we 



find at 24.7 meV. Furthermore, as we commented in the 
previous section, the 11.3 meV TA2 resonance has atomic 
displacements parallel to the surface plane, and would 
normally be undetectable by neutral-atom scattering. 

These authors also report a study of the surface dy- 
namics using a central nearest-neighbor harmonic force 
constant model. Qualitatively, this simple model cor- 
rectly reproduces the zone-center surface excitations 
which we have calculated. The agreement between the 
two theories is at first surprising in view of the rather 
large changes of the surface harmonic interactions which 
we presented in Tables 2 and 3. However, in the long 
wavelength limit these surface features are "leaky" exci- 
tations which penetrate well below the first atomic layer 
of the surface. The calculated frequency of such an exci- 
tation is relatively insensitive to the precise value of the 
surface interatomic interactions, and instead depends on 
the fact that propagating elastic waves will be perfectly 
reflected from the surface. Close to the surface Brillouin 
zone boundaries the effective penetration depth of these 
modes is substantially smaller, and the frequencies are 
much more sensitive to the surface interatomic poten- 
tials. 

More recently Durukanoglu, Kara and Rahman (DKR) 
reported a theoretical study of the equilibrium struc- 
ture and dynamics of Cu vicinals using the embedded 
atom method (EAM) For Cu(211), these authors 

also observed enhanced structural relaxations near the 
step edges as we do, and the magnitudes of their pre- 
dicted relaxations are consistent with the trends we dis- 
play in Table I. They also compute interatomic force 
constants near the step edges, and find that the inter- 
atomic force constants are enhanced near the step edges. 
DKR also note that the diagonal elements of the dynam- 
ical matrix are decreased relative to their bulk values; 
this follows directly from the lower coordination of the 
surface atoms. Thus the off-diagonal terms provide a 
useful quantitative measure of the enhancement of the 
"back bonding" of under-coordinated atoms near the step 
edges. Although the EAM and LDA results both indicate 
an enhancement of these interactions, the enhancement is 
approximately twice as large within the LDA than that 
obtained with the EAM model. The local density ap- 
proximation provides a more complete description of the 
relaxation of the conduction charge density, particularly 
near the step edges. In our view this quantitative dis- 
crepancy between the calculated force constants reflects 
the fact that these depend rather sensitively on the re- 
laxation of the charge density near the step sites. 

Our calculations point to an important piece of new 
physics at the stepped surface. As we emphasized in 
Section VI, the surface lattice dynamics is dominated 
by a simple zone folding effect by which the flnite mo- 
mentum Rayleigh wave, and longitudinal resonance are 
folded back through the surface Brillouin zone center F. 
Normally these excitations are observed as modes with 
flnite momentum transfer in the surface plane and are 
accessible only by neutral atom scattering or by inelastic 
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electron scattering in the impact regime. On the vicinal 
surface, these backfolded modes can be coupled to charge 
fluctuations near the surface, an effect which we expect 
to be particularly strong near the step edges. Thus these 
excitations can be coupled to a macroscopic charge fluc- 
tuation of the surface dipole layer, and are accessible by 
inelastic electron scattering. This is particularly signifi- 
cant because inelastic electron scattering can be used to 
study the dynamics of the conduction charge near step 
edges, and to examine the effects of adsor bates, step wan- 
dering, etc. on this coupling. These are central issues in 
the surface chemistry associated with the step edges, and 
we are therefore currently investigating this effect theo- 
retically by calculating the strength of this coupling. 

In summary, we have reported the first ab-initio study 
of the structure and dynamics of a stepped surface for a 
transition metal. We find that atomic relaxations near 
the step edges for the Cu(211) surface are substantially 
enhanced with respect to their counterparts at a low in- 
dex surface. This produces very large changes of the 
surface harmonic force constants, much larger than in- 
ferred from previous experiment and theory. For the 
Cu(211) surface we find that there are three important 
step-derived excitations in the Brillouin zone center vi- 
brational spectrum. Our result for the frequency of 
one mode, corresponding to a backfolded Rayleigh wave 
agrees well with experimental measurements using inelas- 
tic He scattering. The two othcir features were inacces- 
sible to these measurements, and we have suggested a 
method for observing these other excitations. 

This work was supported by the NSF under grants 
DMR 93 13047 and DMR 97 02514. 
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TABLE II. Selected intcrlaycr surface harmonic interaction 
constants are compared with their values in the bulk. The 
atoms and the polarization of the displacements are labeled 
according to the notation of Figure 1. All constants are in 
units of eV/A^. 





Polarization 


i!'aa(l,4) 


K„c«(l,2) 


surface 


a = X 


4.52 


0.41 


bulk 


a = X 


2.89 


0.12 


surface 


a = y 


1.47 


1.27 


bulk 


a = y 


0.92 


0.92 


surface 


a = z 


-0.49 


3.24 


bulk 


a = z 


-0.21 


2.70 
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TABLE III. Diagonal components of the harmonic force 
constants for the surface terrace atoms on Cu(211). The val- 
ues of Kxx{i, i) are given in eV/A'^ for the three atoms along 
the surface terrace, using the labeling convention of Figure 1. 
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FIG. 2. Atomic displacements in the equilibrium structure 
of our model thin film for Cu(211). The arrows show the 
magnified displacements of each atom from its position in the 
ideal lattice termination to the equilibrium position. For clar- 
ity the displacements have been magnified by roughly a factor 
of 20. 



FIG. 1. Structure of the Cu(211) vicinal surface. The sur- 
face consists of three-atom (111) terraces separated by single 
atom (100) steps. The lattice parameters and conventions for 
our coordinate system are illustrated in the top panel. We 
model the (211) surface with a film containing 17 Cu sites. 
The lower panel is a projection of this film on the (Oil) plane 
with the atoms in the unit cell identified. 
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FIG. 3. Vibrational densities of states for the Cu(211) film. 
The spectra give the densities of states for modes which are 
even under reflection through the (Oil) mirror plane. The 
light curve is a trace over even modes of the expanded film. 
The bold curve is a trace over the modes weighted by the 
squares of the normal-mode amplitudes in the surface terrace. 
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FIG. 5. A snapshot of the displacement field for the back- 
folded Rayleigh mode at 11.5 meV in the surface vibrational 
spectrum. The amplitudes decay exponentially into the bulk. 
At the surface, this mode has a standing wave pattern which 
is "pinched" at the step edges. 



FIG. 4. Vibrational densities of states for the Cu(211) film. 
The spectra give the densities of states for modes which are 
odd under refiection through the (Oil) mirror plane. The light 
curve is a trace over odd modes of the expanded film. The 
bold curve is a trace over the modes weighted by the squares 
of the normal-mode amplitudes in the surface terrace. 
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FIG. 6. A snapshot of the displacement field for the back- 
folded longitudinal resonance centered at 24.7 meV. The mode 
is weakly resonant with bulk excitations. At the surface, the 
atomic displacements are mainly polarized parallel to the sur- 
face. 




FIG. 7. A snapshot of the displacement field for the TA2 
resonance centered at 11.3 meV. This mode is odd under re- 
flection through the (Oil) mirror plane, and can be described 
as a shear horizontal mode polarized parallel to the surface 
steps and localized at the surface step edge. 
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